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Investigation of the Two-Particle-Self-Consistent Theory for the 
Single-Impurity Anderson Model and an Extension to the Case of 

Strong Correlation 

Tetsuro SASO0 

Department of Physics, Faculty of Sciences, Saitama University, Urawa 338-8570 
(Received ) 

The two-particle-self-consistent theory is applied to the single-impurity Anderson model. It is 
found that it cannot reproduce the small energy scale in the strong correlation limit. A modified 
scheme to overcome this difficulty is proposed by introducing an appropriate vertex correction 
explicitly. Using the same vertex correction, the self-energy is investigated, and it is found that 
under certain assumptions it reproduces the result of the modified perturbation theory which 
interpolates the weak and the strong correlation limits. 
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§1. Introduction 

The single impurity Anderson model (SIAM)S is not 
only a modeL for the magnetic impurity in metals, 
chemjsorptiontf and a quantum dot in mesoscopic sys- 
tems,™ but also the most important model for the study 
of the strongly correlated electron systems. It was proved 
that the lattice models with strong correlation is reduced 
to a generalized impurity Anderson model with energy- 
dependent hybridization in the limit of infinite spatial 
dimensions .[§ In that theory, a knowledge on the excita- 
tion spectrum of the impurity is indispensable. 

It is, however, not possible to obtain an exact solu- 
tion to the dynamical properties of SIAM. The best way 
known to date is to-jresort to the numerical renormaliza- 
tion group method,© but it needs heavy numerical tasks, 
especially when both finite temperature, finite magnetic 
field and finite excitation energy must be taken into ac- 
count. Therefore, it is desirable to find an approximate 
but a simpler scheme which is applicable to the problems 
in realistic materials with strong correlation. 

Various approximate schemes proposed to date have 
several deficiencies despite theni advantages. The non- 
crossing approximation (NCA)ErLrEf cannot reproduce 
the Fermi liquid properties at low frequencies and low 
temperatures, and the slave boson mean field or l/N 
expansion theoryEP cannot avoid redundant phase tran- 
sitions. The second-order perturbation theory exhibits 
Fermi liquid properties and happens to reproduce the 
correct atomic limit in the symmetric casejiif but it is 
not the case when there is an electron-hole asymmetry. 
This point is recently overcome by introducing an inter- 
polative form for the self-energy phenomenologically.El 

Therefore, it is of much importance to establish a 
method to calculate the static and dynamic properties of 
SIAM in the wide range of parameters efficiently. Such 
a method must be free from the Hartree-Fock instabil- 



ity and have proper energy scales. In the present pa- 
per, we first-apply the two-particle self-consistent the- 
ory (TPSC)EJ to SIAM in §2. TPSC has proved to 
be successful in the study of the Hubbard model in the 
weak and the intermediate correlations, and is consid- 
ered to be superkir to the fluctuation-exchange approx- 
imation (FLEX)o' in that only TPSC can reproduce 
the side peaks in the single-particle spectrum— Appli- 
cation of FLEX to SIAM is already reported^ which 
has clarified that FLEX is poor on magnelki susceptibil- 
ity than the simple perturbation theory,ES in addition 
to the above-mentioned deficiency on the spectrum. De- 
spite the superiority of TPSC, we will find below that 
it cannot reproduce a small energy scale in SIAM in the 
strong correlation limit. 

TPSC can be also regarded as an extension of 
the selffCpnsistently renormalized spin-fluctuation theory 
(SCR) .113 The-author previously proposed an extension 
of SCR theory^ to-the case of the strong correlation on 
a microscopic basis S3 In §3 we apply a similar scheme 
to SIAM in combination with TPSC instead of SCR, and 
establish a scheme which can be applied to the case of 
strong correlation. We also investigate the self-energy 
using the same vertex corrections and find that under 
certain assumptions our theory -cepxoduces the modified 
perturbation theory (MPT) JilB Ef which interpolates 
the weak and strong correlation limits. Thus the phe- 
nomenological character of MPT is partially resolved. 

§2. The Two-Particle-Self-Consistent Theory 
for the Single-Impurity Anderson Model 

The Hamiltonian for SIAM is written as 

H = ^ 6 k c ta C kcr + E d ^ n d<y 
fccr iy 

+V^2(d+c krT +c+ ij d a ) + Un dl n dl (2.1) 
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in the ordinary notation, where ek and denote the 
energy of the conduction and d electrons, respectively, V 
the hybridization, and U the Coulomb repulsion between 
d-electrons. In the following, the subscript d in rida will 
be dropped. 

In TPSC, the charge and the spin susceptibilities are 
expressed in the form similar to those of RPA, but with 
the renormalized interaction parameter U c and U s for 
charge and spin channels, respectively. These parameters 
are determined by the sum rules, 

2IIoM 



T 



^ 1 



2n (kj) 

- U s Uo(iuj) 



= n + 2(n T n x ) -n , (2.2) 



= n- 2(n T n x ). 



(2.3) 



Here ui denotes the Matsubara frequency, n = (n^) + (nj,) 
and U s are related to {n-^n^) by U s — U{n^n{) / (n-|-)(nj,). 
Thus, U s can be determined from the second equation, 
while U c can be obtained by solving the first equa- 
tion once U s is solved. We assume the symmetric case 
Ed = —U/2 and the paramagnetic state (n a ) = n/2 
throughout the present paper. An extension to the asym- 
metric case is straightforward. Ho(lo) denotes the polar- 
ization function for U = 0, which is calculated at T = 



noM 



2A 



log (1-5 



(2.4) 



ttu>(uj + 2iA) 

where A denotes the resonance width of the impurity, 
A = irp c V 2 (p c is the density of states of conduction 
electrons at the Fermi energy). At low frequencies, it 
behaves as 

rioH ~ n ° (u> « A), (2.5) 

1 — iu>/ A 

1/ttA. Using this form, the sum rule eq.(2.3) 



with Ilo 
reads 



U s 2 uj c 



(2.6) 



where oj c denotes a cutoff frequency of the order O(A) 
and T = A - U s /w. 

For U — » 00, we have (n^ni) — ► and hence 
U s — > 0. Then the above equation is solved as T = 
u> c exp(— 7r 2 /2) « w c /139. Namely, we obtain 



Us = 7r[A-w c e-" 2 / 2 ]. 



(2.7) 



If we set gp,B '- 
susceptibility \ 



1 (g = 2 is the g-factor), the static 
X zz is given by 



X(w = 0) = - 



1 



2 1 



duced at (7 > A so as to avoid the magnetic insta- 
bility, but does not reach irA, namely, U s — > 7r[A — 
w c exp(— 7r 2 /2)] and hence the Kondo temperature re- 
mains finite even at f = 00. Therefore, TPSC cannot 
describe the case of strong correlation properly. We will 
seek for an improvement of TPSC in the next section. 
Before it, however, we briefly discuss on the self-energy 
in TPSC. 

The self-energy for the d-electrons is calculated asE3 



£ ff (ie) = £/n_ CT + jT^ 



2[/ s n (iw) 



2t/ c n M 



1 -C/JIoM 

i + . c n„Mj G °( i£ + i ^ < 2 - 9 > 

In the weak coupling limit C/ — > 0, this is reduced to the 
second-order self-energy 

E„.(ie) = Un- a + U 2 T^2u (iLj)G {ie + iw). (2.10) 

It is well known, however, that the second-order self- 
energy does not have an expected spectral properties in 
the electron-hole asymmetric cases. For example, the po- 
sitions of the Ed and Ed + U peaks are not correct, and 
the Friedel sum rule is not satisfied. These difficulties 
can be most easily overcomeJpy the use of the interpola- 
tive form of the self-energyM* and by-subtracting £(0) 
from £(e) as pointed out by Yamada.Ej) More improved 
treatmeuLis-possible by introducing the effective d-level 
positioniM 

§3. Improvement of Two-Particle-Self-Consistent 
Theory 

The dynamical susceptibility of the Anderson impurity 
can be generally expressed by the diagram shown in Fig.l 
and the equation, 



X 



+- 



(iw) = T^x + (ie, ie', i^), 



(3.1) 



X + -(ie,ie',i(j) = -G T (fe)G ; (ie' + iw) 

x 6 e , € , + TJ2 r(ie", ie', iu)x + ~(ie", ie' , & 



(3.2) 



We calculate the irreducible vertex function r(ie, ie', ito) ap- 



= . (2.8) 

C/ S n 27rr v ; 

Therefore, if we define the Kondo temperature Tj^ by 
x (uo = 0) = 1/4T K , we obtain T K = {%/2)T = 
(tt/2)ui c exp(— 7r 2 /2), which is finite even for U — > 
00. Of course, the., correct formula is Tj^ = 
v /[/A/2exp(-7rt//8A)j£3 ) which vanishes for U -> 00. 
For U < 00, a similar analysis yields Tj^(U) = 
(7r/2)cj c cxp[-(vr 2 /2)(l - U S /2U)], which is larger by 
exp[(7r 2 /2)(t/ s /2C/)] than T K (U = 00). 

Thus, in TPSC, the renormalized interaction U s is re- 




e ' +co e+co e+co e"+co e'+co 




Fig. 1. The Feynman diagram for the susceptibility \ + (ie, ie' , ioj) . 
The square denotes the vertex function r(ie, ie', iu>). 



proximately by taking account of the maximally crossed 
diagrams (Fig. 2) as 



U 



r(ie,ie',ia;)= 1 + ^ ( , £+ , e 



(3.3) 
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Fig. 2. The diagram for the vertex function r(ie, ie', iui). 



K(ie + ie',iuj) = T^Gx(ie + ie'-ie")G T (ie" + iw). (3.4) 

e" 

We use the unperturbed Green's function Go a for G CT and 
denote K by Kq. Then, we find that Ko(Q,u>) — IIo(u/), 
Kq(oj, 0) = IIq(cj) (note that e + e' is bosonic) and thus 
at low frequencies 



K (0,u) = K (u,0) 



1 1 
7rA 1 - iu/A' 



(3.5) 



If we approximate r(ie, ie', iuj) by r(io>) = T(0, 0, iw), then 
we obtain 



n (iw) 



i-r(iw)n (iw)' 



r(iw) = - 
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(3.6) 



(3.7) 



cm (iw) 

Note that at T = and the low frequency limit, the 
effective interaction becomes 

U 



r(o) = 



(3.8) 



l + U/irA' 

which smoothly interpolates between T(0) = U for small 
U and r(0) — » 7rA for [7 — > oo. Therefore, the magnetic 
susceptibility x Z2 (0) = (l/2)x H (0) diverges when and 
only when U — > oo. This is a desired property. Within 
the present approximation, eq.( |3.6D becomes 

x + -(iu) = noM[i + cmo(iw)]. (3.9) 

Namely, the terms higher than the second order with re- 
spect to U in the perturbational expansion of x( w ) are 
completely canceled out with the terms from the ver- 
tex correction, which is of course an artifact of the ap- 
proximation. Note also that eq. (|3.9| ) is correct up to 

o(u)Em 

Using the low frequency form for IIo(u;), we obtain 



x + ~M 



2xo 



l-iw/A' 



(3.10) 



where xo = (l/27rA)(7rA + U)/tvA and A = A/[l + 
U/(irA + [/)]. On increasing U from to oo, Xo varies 
from l/(2vrA) to U/2(nA) 2 -> oo. Thus the Kondo tem- 
perature Tj^ — 1/ 4xo vanishes for U — > oo as we desired, 
but more slowly as — > (irA) 2 /2U than the correct be- 
havior. This must be improved in the future by taking 
account of better vertex corrections. On the contrary, A 
is reduced from A only to A/2 for U = — > oo, and re- 
mains finite at U = oo. Therefore, we replace A with an 
effective parameter A' and determine it so as to satisfy 
the sum rule eq.(2.3) with U s = T(0). Namely, in the 



symmetric case and at T = 

1 r(o) = i 

2U 7T 



dtu Im 



4xo 



10 

4xoA' 



1-iu/A' 



log 



(3.11) 



For U — 0, we have T(0) = 0, and the cutoff lo c should 
be chosen as lo c = Aexp(7r 2 /4) since xo — 1/(2ttA) and 
A' = A. For U > 0, the cutoff may also be chosen 
as a similar form ui c = A' exp(7r 2 /4). Then we obtain 
A' = A x 27rA/(7rA + U), which vanishes as U increases. 
If we use eq.(B^) for H (uj), the integration converges 
without a cutoff. In that case, we previously proposedH 
to modify ( CJ ) hito 



X 



1 



(3.12) 



and determine the parameter C by the sum rule. This is 
essentially corresponding to modifying A as mentioned 
above. At low frequencies, x" 1 behaves as 

2xo 



X 



where 



A' = (A- 



l-iw/A'' 



2Gxo) 



(3.13) 



(3.14) 



which tends to A(7r/C)7rA/(7rA + U) -> for U -» oo. 
In Fig. 3, we show this A' and compare it with A' = 
A x 2ttA/(ttA + U) in the above analysis using eq.(3.11). 
Both are in good agreement for U 3> A and become 
smaller than A for U/A>7 although the reduction of 
A' is not sufficient quantitatively. The first rise of A' at 
U^A is to compensate the decrease of the left-hand side 
of the sum rule 



1 



duj Im- 



2IT0M 



r(o) 



= 1 



(3.15) 



1 - T(lu)U (ui) 2U 

from 1 to 0.9 for U/A — to 2 when the C term is 
not included (see Fig. 4). This decrease is traced back 
to the decrease of r(0)/2C7 = (l/2)/(l + U/ttA). For 
U/A > 2, the increase of the first term of the l.h.s. of 
the above equation dominates and l.h.s. increases again. 
The imaginary parts of the dynamical susceptibility cal- 
culated with this C term correction are shown in Fig. 5. 
We set A = 1 here and henceforth. 

Next, we investigate the self-energy. The Feynman 
diagram is shown in Fig. 6, which may be expressed as 

E T (ie) =U ni + --- 

+ T 3 ]T r(ie ! ie , ,i W )G T (ie , )G;(ie , +ia;) 

x r(ie', ie",iw) • • • G T (ie")Gx(ie" + iu) 

x r(ie",ie,iw) xG;(ie + iw) + ---. (3.16) 

Keeping the behaviors of K o (0,oj) and K o (e,0) in mind, 
we approximate r(ie, ie', iuj) at the left and the right 
ends by the separable form [f '(ie)r(iw)] 1 / 2 , where f'(ie) 
denotes some average of r(ie,ie',iw) on e' and uj. To 
fix f'(ie), we consider the first order self-energy due to 
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Fig. 3. A'/A calculated by eq.(3.14) is shown as a function of 
U. Also shown are A'/A = 2wA/(nA + U) (dashed line) and 
A/A = 1/[1 + U/{nA + U)] (dotted line). 



Fig. 6. A typical Feynman diagram for the self-energy. 



1.5 



0.5 



4 6 
U/A 
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Fig. 4. The left-hand side of eq.(3.15) is plotted by the full line 



for the case without the C term correction. The first term in 
l.h.s. is plotted by the dashed line, whereas the second term by 
the dash-dotted line. 




CO 



Fig. 5. The imaginary parts of the dynamical susceptibility 
X + ~(^Y are shown for U/A =0, 2, 4 and 8. 



T(ie,ie',iw) (Fig.7): 

S T (ie) = T]T r(ie, ie, ico)G i (ie + iw), (3.17) 

which is expanded for small U as 

= TJ2[U- U 2 K (ie + fe.iw) + ■ ■ •] G x (ie + iu). 



f/nx + E; 2) (ie) + ---. 



(3.18) 



where Si (ie) denotes the second-order self-energy. On 
the other hand, if we replace r(ie,ie,iw) in Sf(ie) with 
the average f (ie) = (r(ie, ie, 



S T (ie) = f(ie)T^G oi (ie + i^) 

= Un l - U 2 n l (K (ie + ie, iw)) w + 

Comparing these two expressions, we obtain 

(Ko(ie + ie,iLj)} u = -•E<f>{i £ )/U 2 n l . 
Therefore, we approximate f(ie) as 

f(ie) 



U 



l-B^Y>(ie)/Uni 



(3.19) 

(3.20) 
(3.21) 

(3.22) 



where B' is an adjustable parameter to take account of 
the effects of the higher order terms. Approximating 
f '(ie) by T(ie) and replacing the intermediate T(ie, ie', iuj) 
by the previous T(iw), we obtain 



E T (ie)~[/n i + 



f(ie)T£l» T 



n (iw) 



- r(iw)iioM 



Goi(ie + iuj). (3.23) 



We further note that 



r(ic) 



u, 



(3.24) 



i - r(iw)n (iw) 

which holds within the present approximation. Therefore 



Two-Particle Self-Consistent Theory for the Single-Impurity Anderson Model 



5 



we obtain 

E T (ie) 



Urn + 



Urn 



^ 2 rE w n MGoi(ie + i^) 



B'Hf>(ie)/Ui H 



sf>(ie) 



l-SE^ie)' 



(3.25) 



where B = B' /Uni- This is exactly the same form as the 
inteipolative self-energy proposed by Martin-Rodero, et 
al.,El which bridges between the second-order and the 
atomic limit self-energies in both the symmetric and the 
asymmetric cases. B was determined so as to reproduce 
the correct self-energy in the atomic limit when U — > oo. 
Namely, one obtains B = (1 — 2?i_ (T )/[/?i_ (J (l — n^ a ). 




Fig. 7. The left- and right-hand sides denote eqs.(3.19) and 
(3.18), respectively. 



Since the above formula for the self-energy does not 
satisfy the Friedel sum rule, one has to introduce an ef- 
fective jL level energy to be adjusted to -fulfil the sum 
rulejSEf or to subtract S(0) from S(e)t£l ) to approx- 
imately satisfy it. Details of the calculations will be 
found in refj2^. Note that B vanishes for the symmet- 
ric case. Namely, the self-energy becomes equal to the 
second-order one due to the cancelation of the higher or- 
der terms in the present approximation. Recently, it was 
pointed out that the second-order self-energy becomes 
exact-when the width of the conduction band becomes 
zero.eS The cancelation of the higher order terms in our 
treatment may be related to this finding. 

§4. Conclusions 

In this paper, we have applied the two-particle self- 
consistent approximation to the single-impurity Ander- 
son model, and found that it is not suitable for the de- 
scription of the strong correlation limit. The form of the 
dynamical susceptibility was assumed to be of RPA type 
in TPSC, and the effective interaction U s is determined 
by the sum rule eq.(2.3). In this scheme, a magnetic in- 
stability was avoided but U s remained to be too small in 
U — > oo. Thus, the low energy scale remains finite and 
does not vanish even at this limit. This is understandable 
since TPSC was devised originally for the cases of weak 
and intermediate correlations. We needed an improved 
treatment of the vertex corrections for the dynamical 
susceptibility. 

Then, we have proposed a modification of TPSC appli- 
cable to the strong correlation. We have used the renor- 
malized interaction T(uj) = U/[l + UHq(lo)] instead of 



U s , and assumed the modified form of x( w ) as 

X +-( W ) = ™M , (4.1) 

x K ' i-rHnoH-2iCwrioM' v ; 

determining the parameter C by the same sum rule. This 
scheme has realized the vanishment of the energy scale 
in the strong correlation limit, although the decrease of 
the energy scale is slow compared to the correct behavior 
(Kondo temperature). Apparently, this is due to our use 
of the approximate vertex corrections to x( w ) which is 
not exactly consistent with the self-energy diagram. It 
is known that the calculation of the static susceptibility 
by the numerical differentiation of the electron numbers 
n-f — ni within MPT leads to a rather^prrect values if it 
is combined with the Friedel sum ruleH But one cannot 
apply it to the calculation of the dynamics. 

We have also investigated the self-energy within the 
same vertex correction and found that with an appro- 
priate approximation it reproduces the so-called modi- 
fied perturbation theory which interpolates the perturba- 
tive regime and the atomic limit. Although the present 
derivation includes a crude approximation, it suggests 
that an appropriate treatment of the vertex correction 
may lead to a theory which bridges the weak and strong 
correlation limits. Extension of the present scheme to 
the asymmetric case is straightforward. In order to im- 
prove the present theory further, we need a better vertex 
function. It would be of much interest to calculate the 
set of eqs.(3.1-4) exactly. 

We have already applied a theory, which is similar to 
the present one but is combined with SCR theory, Jtp 
the description of the quantum critical phenomena. i3 
It will be improved by the application of the present 
formulation to lattice problems. Such a study is now 
in progress, and will be useful for the investigation of 
the realistic materials with strong correlation when it is 
combined with the band calculations. 
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